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ABSTRACT 

We present the magnetic and electric field structures Ets well as the currents and charge 
densities of pulsar magnetospheres which do not obey the ideal condition, E ■ B — CL Since the 
acceleration of particles anti the production of radiation requires the presence of an electric field 
component parallel to the magnetic field. Eg, the structure of non-ideal pulsar magnetospheres 
is intimately related to the production of pulsar radiation. Therefore, knowledge of the structure 
of non- Ideal pulsar magnetospheres is important because their comparison (including models 
for the production of radiation) with observations will delineate the physics and the parameters 
underlying the pulsar radiation problem. We implement a variety of prescriptions that support 
nonzero values for Ey and explore their effects on the structure of the resulting magnetospheres. 
We produce families of solutions that span the entire range between the vacuum and the (ideal) 
Force- Free Electrodynamic solutions. We also compute the amount of dissipation as a fraction 
of the Poynting flux for pulsars of different angles between the rotation and magnetic axes and 
conclude that this is at most 20-40% (depending on the non-ideal prescription) in the aligned 
rotator and 10% in the perpendicular one. We present also the limiting solutions with the 
property J ^ pc and discuss their possible implication on the determination of the “on /off” 
states of the intermittent pulsars. Finally, we find that solutions with values of J greater than 
those needed to null Ey locally produce oscillations, potentially observable in the data. 

Subject headings: Magnetohytlrodynamics (MHD)-Methods: mi merical- pulsars: general- Gamma rays: 
stars 


1. Introduction 

Pulsars are extraordinary objects powered by 
magnetic fields of order 10 12 G anchored onto neu- 
tron stars rotating with periods ^ 10” 3 — 10 s. 
These fields mediate the conversion of their rota- 
tional energy into MHD winds and the accelera- 
tion of particles to energies sufficiently high to pro- 
duce GeV photons. Their electromagnetic emis- 
sion is quite complex and ranges from the radio to 
the multi-GeV 7— ray regime. 

Pulsar radiation has been the subject of many 
studies since their discovery; however, despite 
such a multi-decade effort by a large number of 
researchers its details are not fully understood. 


There are several reasons for this. On the ob- 
servational side, the sensitivity of the 7 — ray in- 
struments (the band where most of the pulsar ra- 
diation power is emitted) has become adequate 
for the determination of a significant number of 
pulsar spectra only recently. Since its launch in 
2008, the Fermi Gamma-Ray Space Telescope lias 
revolutionized pulsar physics by detecting nearly 
90 7— ray pulsars so Far, with a variety of well- 
measured light curves. On the theoretical side, 
the structure of the pulsar magnetosphere, even 
for the simplest axisy mine trie case, remained un- 
certain for a long rime. The modeling of pulsar 
magnetospheres lias seen rapid advancements only 
very recently. 


t 



The original treatment of the non -axisymmetric 
rotating stellar dipole in vacuum (Deutsch 1955) 
established the presence of an electric field com- 
ponent, B||\ parallel to the magnetic field on the 
stars surface. Then Goldreich & Julian (1969) 
suggested that tins component of the electric field 
would pull charges out of the pulsar, opening up 
magnetic field lines that cross the Light Cylinder 
(hereafter LC), to produce an MHD wind with cur- 
rents flowing out of and into the neutron star polar 
cap. The charge density needed to achieve this is 
known as the Goldreich - Julian density />gj ; an as- 
sociated current density Jqj = cpc j is obtained if 
we allow these charges to move outward at nearly 
the speed of light* 

Following the work of Goldreich k Julian 
(1969), a number of attempts were made to pro- 
duce the structure of the pulsar magnetosphere. 
Thus, Scharlemann k Wagoner (3973) were able 
to reduce the structure of the axisymmetric mag- 
netosphere to a single equation for the poloidal 
magnetic flux, the so-called pulsar equation. How- 
ever an exact solution of tills equation that would 
produce the magnetic field structure across the 
LC was missing for almost 3 decades (for a review 
of the problems and the various attempts to pro- 
duce a consistent pulsar magnetosphere, even in 
the axisymmetric case* see Michel 1982). In fact, 
the inability of models to determine the magne- 
tos p her ic structure across the LC led people to 
speculate that its structure is not smooth across 
tills surface and assumed that the observed radia- 
tion is the result of magnetic field discontinuities 
across the LC (Mestel k Shibata 1994). 

The structure of the axisymmetric pulsar mag- 
netosphere within Force Free Electrodynamics 
(hereafter FFE) (he* assuming that the only 
surviving component of the electric field is that 
perpendicular to the magnetic field, or equiva- 
lently that E B 0) was given by Contoponlos, 
Kazanas* k Fendt (1999) (hereafter CKF) who 
used an iterative procedure to determine the cur- 
rent distribution on the LC so that the magnetic 
field lines should cross this surface smoothly. Be- 
sides allowing a smooth transition of the magnetic 
field lines across the LC. this solution also resolved 
the issue of the current closure in pulsar magne- 


Ei; = K B)B/B 2 , where E, B an? the electric and mag- 
netic field vectors. 


tospheres: The current starts near the polar axis 
and closes, for ft > /?lc (here /£lc is the radius of 
the LC), mainly along an equatorial current sheet 
(allowing also a small amount of return current at 
a finite height above it). For R < Rt c the Row is 
mainly along the separatrix of the open and closed 
field lines* The axisymmetric solution has since 
been confirmed and further studied by several oth- 
ers (Gruzinov 2005; Timokhin 2006; Komissarov 
2006; McKinney 2006), in particular with respect 
to its properties near the Y-point (the intersection 
of the last closed field line and the LC; Uzdensky 
2003). 

The first non-axisymmetric (3D), oblique rota- 
tor magnetosphere was presented by Spitkovsky 
(2006), who used a time-dependent numerical code 
to advance the magnetic and electric fields under 
FFE conditions to steady state. These simulations 
confirmed the general picture of current closure 
established by the CKF solution and produced a 
structure very similar to that of CKF in the ax- 
isymmetric case. Similar simulations were per- 
formed by Kalapotharakos k Coutopoulos (2009) 
who, using a Perfectly Matched Layer (hereafter 
PML) at the outer boundary of their computa- 
tional domain, were able to follow these simula- 
tions over many stellar rotation periods. In gen- 
eral, the 3D magnetosphere, just like the axisym- 
metric one, consists of regions of closed and open 
field lines with a large scale electric current circuit 
established dong open magnetic field lines. In the 
3D case, the current sheet needed for the global 
current closure is in fact undulating, as foreseen 
in the kinematic solution of Bogovalov (1999). As 
shown by recent simulations* the undulating cur- 
rent sheet structure is stable at least to distances 
of 10/? and for several stellar rotation periods 
(Kalapotharakos* Coutopoulos. k Kazanas 2011). 

The FFE solutions presented so far, by con- 
struction. do not allow any electric field compo- 
nent parallel to the magnetic field B E B 0): 
as such they do not provide the possibility of 
particle acceleration and the production of radia- 
tion* in disagreement with observation* Nonethe- 
less. models of pulsar radiation have been con- 
structed assuming either a vacuum dipole geome- 
try, or the field line geometry of the FFE magne- 
tosphere. Polar cap particle acceleration, - -ray 
emission and the formation of pair cascades were 
discussed over many years as a means of account* 
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mg for the observed high-energy radiation (Rud- 
erinan k Sutherland 1975: Daugherty k Harding 
1996). However, measurements by Fermi of the 
cutoff si jape of the Vela pulsar spectrum (Abdo 
et ah 2009b) has ruled out the super-exponential 
shape of the spectrum predicted by polar cap 
models due to the attenuation of pair production. 
Models placing the origin of the high-energy emis- 
sion in the outer magnetosphere or beyond are 
now universally favored. Slot gap models (Mus- 
limov k Harding 2004; Harding et ah 2008), with 
acceleration and emission in a narrow layer along 
the last open field line from the neutron star sur- 
face to near the LC, and outer gap models (Ro- 
mani k Yadigaroglu 1995; Romani 1996; Takata 
et ah 2007) produce light curves and spectra in 
first-order agreement with the Fermi data (Ro- 
mani k Watters 2010; Venter, Harding, k Guille- 
mot 2009). All of the above models assume a 
retarded vacuum magnetic field. There are also 
kinematic models of pulsar emission derived by 
injecting photons along the magnetic field lines of 
either the FFE magnetosphere (Bai k Spitkovsky 
2010b) or those of the vacuum field geometry (Bai 
k Spitkovsky 2010a), or along the electric current 
near the equatorial current sheet where the ratio 
Jj pc — > 1 in FFE (Contopoulos k Kalapotharakos 
2010). These models also produce light curves In 
broad agreement with observations. However, fits 
of Fermi data light curves with slot gap and outer 
gap models in FFE geometry are less favorable 
than those for the vacuum dipole geometry (Hard- 
ing et al. 2011). 

While FFE magnetospheres do not allow for 
particle acceleration, they do, however, provide 
a global magnet ospheric structure consistent with 
the boundary conditions on the neutron star sur- 
face (a perfect conductor), a smooth transition 
through the LC and the establishment of a large 
scale MHD wind. They also determine both the 
currents and the sign of the charges that flow in 
the magnetosphere* However, they do not pro- 
vide any information about the particle produc- 
tion that might be necessary ro support the under- 
lying charge conservation and electric current con- 
tinuity. This was apparent already in the CKF so- 
lution and was also noted in Michel (1982): There 
are regions, most notably the intersections of the 
null charge surface with open magnetic field lines 
and the region near the neutral (Y) point of the LC 


where the FFE charge density changes sign along 
an open field line (along which particles are sup- 
posed to flow freely). The numerical FFE solution 
has no problem in these regions because it can sim- 
ply supply at will the charges necessary for current 
closure and for consistency with the global bound- 
ary conditions. This, however, requires sources of 
electron- positron pairs at the proper regions and 
of the proper strength. Even though the FFE 
treatment of the pulsar magnetosphere involves no 
such processes, these must be taken into account 
sdf-consistentlyv Viewed tins way, pulsar radia- 
tion is the outcome of the “tension” between the 
boundary conditions imposed by the global mag- 
netospheric structure and the particle production 
necessary for charge conservation and electric cur- 
rent continuity. 

In Figs. la,b,c we present the structure of the 
FFE magnetosphere on the poloidal plane - S"2 
for inclination angles a (0°, 45 c ,90°). These fig- 
ures depict the poloidal magnetic field lines and in 
color the value of Jj per The vertical dashed lines 
denote the position of the LC. In the axisym met- 
ric case the value of Jj pc is less than one in most 
space and exceeds this value only in the null sur- 
faces and the separatrix between open and closed 
field lines. The value l of this ratio is significant 
because it denotes the maximum current density 
for carriers of single sign charge. Values of this 
ratio greater than 1, by necessity denote the pres- 
ence of both positive and negative charges moving 
in opposite directions. For a ~ 0 T J/pc < 1 in most 
of the polar cap region, however, as the value of 
inclination a increases so does the size of regions 
with J/pc > 1; in the perpendicular rotator we 
see clearly that Jj pc> 1 over the entire polar cap 
region. 

With the above considerations in mind, the 
question that immediately arises is how would 
the structure of the pulsar magnetosphere and 
the consequent production of radiation be modi- 
fied, If the ideal condition, namely E B = 0 , were 
dropped* This is the subject of t he present work. 
A successful magnetospherie configuration (solu- 
tion) should be able to address the following ques- 
tions: a) Where does particle acceleration take 
place? b) How close are these processes to a 


■ Hereafter, wherever we use the ratio J/pc we consider Its 
absolute value. 


3 


steady state? c) Is the ensuing particle acceler- 
ation and radiation emission consistent with the 
Fermi pulsar observations? 

Unfortunately, while the ideal condition is 
unique, its modifications are not. The evolution 
equations require a relation between the current 
density J and the fields, be* a form of Ohm’s taw, 
a relation which is currently lacking. Meier (2004) 
proposed a generalization of Ohm’s law in plasmas 
taking into account many effects (e.g* pressure, 
inertia, Hall). However, in the present case we 
are in need of a macroscopic expression for Ohm’s 
law which would be also everywhere compatible 
with the microphysics self-consistently, a proce- 
dure that might involve Particle in Cell (PIC) 
techniques, well beyond the scope of this paper. 
For this reason, the present work has a largely ex- 
ploratory character* A number of simple prescrip- 
tions are invoked which lead to magnetospheric 
structures with values of Bj| ^ 0. The global elec- 
tric current and electric field structures are then 
examined with emphasis on the modifications of 
the corresponding FFE configuration. We are 
presently interested only in the broader aspects 
of this problem, deferring more detailed calcula- 
tions to future work* In particular, the general 
compatibility of the E<] magnitude with the val- 
ues needed to produce (through pair production) 
self-consistent charge and current densities from 
the plasma microphysics will not be discussed in 
the present paper. The contents of the present 
paper are as follows: In §2 we discuss the differ- 
ent modifications in the expression for the electric 
current that we apply in order to obtain non-ideal 
solutions* In §3 we present our results and the 
properties of the new solutions. In §4 we give 
particular emphasis to specialized solutions that 
obey certain global conditions. Finally, in §5 we 
present our conclusions* 


2- Non-Ideal Prescriptions 


Its the FFE description of pulsal magneto* 
spheres Spitkovsky (2006) anil Kalapotharakos k 
Coiitopoulos (2009) solved numerically the time 
dependent Maxwell equations 


= — cV x E 
dt 

oe 

= r V x B - 4ffJ 
at 


under ideal force-free conditions 


E B * 0, pE + - JxB = 0, 
c 

where p — V ■ E/(4tt). The evolution of these 
equations in time requires in addition an expres- 
sion for the current density J as a function of E 
and B. This is given by 


, ExB c BVxB-EVx E_ 

3 = cp-E5- + 7Z b~2 B C3) 


B 2 


4tt 


(Gruzinov 1999)* The second term in Eq. (3) en- 
sures that the condition E J_ B (the ideal condi- 
tion) is preserved during the time evolution* How- 
ever, this term introduces numerical instabilities 
since it involves spatial derivatives of both fields. 
Thus, Spitkovsky (2006) implemented a scheme 
that evolves the fields considering only the first 
term of Eq. (3), and at the end of each time step 
“kills” any developed electric field component par- 
allel to the magnetic field (E»[)* In the case of 
non-idea! status, the Ejj = 0 condition does not 
apply anymore and electric field components par- 
allel to the magnetic field are allowed* However, 
while in FFE the prescription that determines the 
value of E| is unique (it is equal to zero), there is 
no unique prescription in the non-ideal case. Be- 
low we enumerate the non- Ideal prescriptions that 
we consider in the present study: 

(A) The above implementation of the ideal con- 
dition, hints at an easy generalization that leads 
to non-ideal solutions: One can evolve Eqs. (1) &: 
(2)* using only the first term of the FFE current 
density (Eq. 3), and at each time step keep only a 
certain fraction b of the E,| developed during this 
time instead of forcing it to be zero. In general 
the portion b of the remaining Eij can be either 
the same everywhere or variable (locally) depend- 
ing on some other quantity (e.g. p* J). As 6 goes 
from 0 to 1 the corresponding solut ion goes from 
FFE to vacuum. In this case an expression for the 
electric current density is not given a priori, and 
J can be obtained indirectly from the expression 


J 


4tt 


(cV x B - 



( 4 ) 


( 1 ) 

( 2 ) 


4 


(B) Another way of controlling E is to introduce 
a finite conductivity a. In this case we replace the 
second term in Eq* (3) by ^tE and the current 


density reads 


( 8 ) 


E x B _ , e . 

J *^-§3“ + (tE H (®) 

Note that Eq. (5) is related to but is not quite 
equivalent to Ohm’s law which is defined in the 
frame of the fluid. Others (Lyutikov 2003; Li, 
Spitkovsky, k Tchekhovskoy 2011) implemented 
a different version closely related to Ohm’s law. 
The problem is that the frame of the fluid is not 
well defined in our problem. Moreover, it is not a 
priori known neither that a single form of Ohm’s 
law is applicable in the entire magnetosphere, nor 
that the conductivity corresponding to each such 
form of Ohm’s law should be constant (e.g, in re- 
gions where pair production takes place). This is 
why we chose a simpler expression, which turns 
out to yield results qualitatively very similar to 
other formulations based on Ohm’s law. As is the 
case in prescription (A) with the parameter b , here 
too, as a ranges from a = 0 to u oo we ob- 
tain a spectrum of solutions from the vacuum to 
the FFE, respectively. We note that even though 
the vacuum and FFE solutions are unambiguously 
defined limiting cases, there is an infinite num- 
ber of paths connecting them according to how 
depends on the local physical parameters of the 
problem (p, £\ B etc.), 

(C) It has been argued (Lyubarsky 1996; Gruzi- 
nov 2007: Lyubarsky 2008) that regions corre- 
sponding to space-like currents (* J/pe. > 1) are 
dissipative due to instabilities related to counter 
streaming charge flows. 3 Moreover, the space- like 
current regions should trace the pair production 
areas (i.e. dissipative areas). Gruzinov (2007) pro- 
posed a covariant formulation for the current den- 
sity that introduces dissipation only in the space* 
likc regions while the time-like regions remain non- 
dissipative. The current density expression in the 
so called Strong Field Electrodynamics (hereafter 
SFE) reads 

T cpE x B + {cV +■ 7 V^) i/2 (BuB + E 0 E) 
J “ B 2 + Eq 

( 6 ) 

where 

B% - El = B 2 - E 2 , B a E a = E ■ B. E a > (J (7) 


2 B 2 + El 

1 D 2 + El 

and a is a function of Eq, B 0 with dimensions of 
conductivity (inverse time). The expression for 
J (Eq. 6) produces either a space-like (ff > 0) 
or a null (a — 0) current. Note that this treat- 
ment precludes time- like currents. Such currents 
do exist hut they are only effectively time- like: 
the current J, along with Ejj, fluctuates contin- 
uously, so that while ./ is locally greater than (or 
equal to) pc its average value is less than that, 
i.e. (J) < pc ♦ This behavior is captured by our 
code and is manifest by the continuously changing 
direction of the parallel component of the electric 
field along and against; the magnetic field direction 
(Gruzinov 2008, 2011), 

In prescriptions (A) and (B) above, care must 
be taken numerically so that the resulting value 
of Ei (defined as |E - E|||) be less than B. Vi- 
olation of tills condition usually happens in the 
current sheet where the value of B goes to zero. 
We note that the presence of the Eo term in the 
denominator of the SFE prescription (Eq. 6) and 
of the prescription of J used by Li et al. (2011) in- 
deed prevents the denominator from going to zero 
and alleviates this issue in the current sheet re- 
gions. However, alleviation of this problem does 
not guarantee that all treatments that include this 
Eq term converge to the same solution or that even 
produce the same dissipation rates in the current 
sheet. Thus, the solutions and in particular the 
dissipation rates of SFE and that of Li et al. (2011) 
can be quite different. 

All the simulations that are presented be- 
low have run in a cubic computational box 
20 times the size of the LC radius Rlc*. he 
-IOElc - ■ * 10/?LCj 3 - Outside this area we imple- 
mented a PML layer that allows us to follow the 
evolution of the magnetosphere for several stellar 
rotations (Kalapotharakos k Contopoulos 2009). 
The adopted grid cell size is 0.02/?lo the stellar 
radius has been considered at r* = 0.3i?ic* and. 
all simulations were run for 4 full stellar rotations 
to ensure that a steady state has been achieved. 


3 Xote that space-ltke currents (J/pc > 1) an* formed exclu- 
sively by counter streaming charge Rows. 




3. Results 


We have run a series of simulations using pre- 
scription (B) that cover a wide range of a values 4 . 
In Fig. 2a we present the values of the Poynting 
flux L integrated over the surface of the star for the 
a = 0° (aligned; red), a = 45° (green) and a = 90° 
(perpendicular: blue) versus a in log- linear scale. 
The two solid horizontal line segments of the same 
color denote in each case (a = 0°, 45° , 90°) the cor- 
responding limiting values of the Poynting flux for 
the FFE and the vacuum solutions, respectively. 
The Poynting flux as a function of the inclina- 
tion angle a for the vacuum and the FFE solutions 
reads 


2 n 2 Sl* . 2 

3 — Sm ° 


Vacuum 


p 2 tl 4 


(9) 


(1 -b sin 2 a) FFE 


where are the magnetic dipole moment and 
the rotational frequency of the star, respectively. 
This means that for a specific value of a the cor* 
responding value of the Poynting flux (spin-down 
rate) of each non-ideal modification should lie be- 
tween the two values of Eq. (9), Thus, Fig. 2a 
shows that for the perpendicular rotator {and any 
case with a ^ 0) the ratio of the Poynting flux over 
that of the vacuum reaches near 1 for low values 
of cr (<r/fl 5 0.01; for a - 0 this is true only for 
a — > 0). In Fig. 2b we plot for these simulations 
the corresponding dissipation power 


JF d = I .l -EdV=-jt I *E\ dV 

rj < r<ra ^<**<**2 

( 10 ) 

taking place within the volume bounded by the 
spheres T\ - QSMlc, i*c just above the star 
surface, and r*a = 2.5/?£,c* versus a, We see 
that the energy losses due to dissipation exhibit 
a maximum that occurs at intermediate values of 
o (<t/U ^ 1) for a = 90° and at higher values of 


■We have found from our simulations that the results of 
prescription (A) art very similar to those of prescription 
(B). \V> therefore restrict the discussion in the rest of this 
section to magnetospheres constructed using prescription 
;B). An explicit reference to results under prescription .. A ) 
is nonetheless made later on in ?hLs section in relation to 
results shown in Fig. 5. 


a as the inclination angle a decreases. The dis- 
sipation power Ed goes to 0 as a goes either to- 
wards 0 or oc since both the vacuum (a = 0) and 
FFE (cr — ► oc) regimes are (by definition) dissipa- 
tionless. We note that despite the fact that the 
Poynting flux (on the stellar surface) varies signif- 
icantly with o - , the corresponding energy loss due 
to dissipation is limited and it never exceeds the 
value ~ 0.14fi 2 fi 4 /c?, The surface Poynting flux 
reflects the spin down rate while the dissipative 
energy loss within the magnetosphere reflects the 
maximum power that can be released as radiation. 
Moreover, when E D is measured as a fraction of 
the corresponding surface Poynting flux it is al- 
ways less than 22% for a = 0°, 15% for a = 45° 
and 8% for a 90°. These maximum fraction val- 
ues occur at lower a values than those that cor- 
respond to the maximum absolute values of Ed 
(Fig. 2b, c). 

We note that all the values presented in Fig. 2 
are the corresponding average values over 1 stellar 
period (4 lh stellar rotation). Moreover, we note 
that the values of E D shown may be a little bit 
smaller than the true ones. This is because our 
simulations introduce some numerical dissipation 
even for the ideal- FFE solutions; for our resolu- 
tion, this produces a decrease in the Poynting flux, 
AL, which at r = 3 can reach 10% of its value 
L on the stellar surface. This means that there is 
some energy loss that can not be measured by the 
J ■ E expression for the values provided by our so- 
lutions. LI ct al, 2011 using higher resolution sim- 
ulations seem to obtain a corresponding decrease 
of ^ 5%, Due to this effect, the Poynting flux 
reduction with r in the dissipative simulations is 
always larger titan its true value. This means that 
the true E& value Is also a tittle higher than com- 
puted. since the energy reservoir is a little larger. 
Nevertheless, the error in Ed is less than A L/L. 
We ran also simulations with half grid cell size and 
simulations with the stellar radius set at 0.45/? ic- 
The Ed values in these cases were at most a few 
percent (mostly for small values of «) of the corre- 
sponding L values lower due to the finer descrip- 
tion of the stellar surface. 

In Fig. 3 we present the magnetospheres struc- 
tures that we obtain on the poloidal plane (ji.f 2) 
using prescription (B) for a ^ 21fl. Each row 
corresponds to the structure at a different incli- 
nation angle a (0 e . 45°.90 c as indicated in the 
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figure. The first column shows the modulus of t lie 
poloidal current J p in color scale together with 
its streamlines. The second and third columns 
show, in color scale, the charge density p and E;| 
respectively, together with the streamlines of the 
poloidal magnetic field B p . We observe that for 
this high a value the global structures are similar 
to those of the FFE solutions (see Contopoulos & 
Kalapotharakos 2010). Even for the a = 0° case 
that has the most different value of the Poynting 
flux from that of the FFE solution, the current 
sheet both along the equator (outside the LC) and 
along the separatrices (inside the LC) does survive 
despite the effects of non-zero resistivity. How- 
ever, we do observe magnetic field lines reconnect- 
ing gradually on the equatorial current sheet be- 
yond the LC. The main differences between these 
solut ions and the FFE ones are highlighted in the 
third column where we plot the parallel electric 
field Ey developed during the evolution of the sim- 
ulation. For a = 0° there is a significant parallel 
electric field component along the separatrices as 
well as over the polar caps 0 . The direction of the 
E|j follows the direction of the current J. Thus, Ej 
within the return current region (the separatrices 
and a small part of the polar cap near them) points 
outward while in the rest of the polar cap region 
points inward. There is also a weaker parallel com- 
ponent Ej | along the equatorial current sheet (riot 
clearly shown in the color scale employed). As 
a increases, the inward pointing E|| component in 
the central part of the polar cap becomes gradually 
offset. During this topological transformation the 
branch of the polar cap area corresponding to the 
return current becomes narrower while the other 
one becomes wider. At a — 90 c , Ey is directed in- 
wards in half of the polar cap and outwards in the 
other half. This topological behavior is similar to 
that of the poloidal current in the FFE solutions. 
For a ^ 0 D there is a clearly visible component of 
E in the closed field lines area as well as along 
the equatorial current sheet outside the LC. 

Figure 4 shows the same plots as Fig. 3 and 
in the same color scale but- for <r ^ 1.5Q. Al- 
though the non-ideal effects have been amplified, 
the global topological structure of the magneto- 
sphere has not changed dramatically. We can still 


°We note that the concepts of the separatrix and the polar 
cap are less! clearly in the case of non- Ideal solutions since 
magnetic field lines close even beyond the LC. 


distinguish traces of the current sheet even in the 
a 0° case, whose Poynting flux is substantially 
reduced from that of the FFE case. As noted 
above, the magnetic field lines close now beyond 
the LC, due to the finite resistivity that allows 
them to slip through the outflowing plasma. The 
parallel electric field component Ey is non- zero 
in the same regions as those seen in the higher 
tr— values; its topology also is not very different 
from that of these cases, but its maximum value 
increases with decreasing cr. We note also that de- 
spite the fact that the dissipative energy losses cor- 
responding to the various values of a are quite dif- 
ferent (see Fig. 2a), the corresponding maximum 
values of Ey are quite similar. 

We have also run many simulations implement- 
ing prescription (A) for various values of the frac- 
tion b of the non -zero parallel electric field com- 
ponent Ey. We have covered this way the entire 
spectrum of solutions from the vacuum to the FFE 
one. The ensemble of solutions that connects these 
two limits seems to be similar to that of prescrip- 
tion (B). It seems that there are pairs of the values 
of the fraction b and the conductivity a that lead 
to very similar (qualitatively) solutions. In Fig. 5 
we plot the poloidal magnetic field lines together 
with the parallel electric field component E|| in 
color scale (similar to the third column of Figs. 3 
and 4) for b ~ 0.75 which has the same Poynting 
flux with the simulation of prescription (B) with 
1 7 = 1.5fl. We see that the magnet os phene struc- 
tures of the two solutions are quite close (compare 
Fig. 5 with the third column of Fig. 4). Also, our 
results are qualitatively very similar to those of Li 
et al. (2011) who, as noted, implemented a much 
more elaborate prescription for J based on Ohm's 
law. 

Figure 6 is similar to Fig. 2 but for prescrip- 
tion (C) - he. SFE, assuming the value of a to be 
constant. This prescription does not tent l to the 
vacuum solution for u 0 but to a special solution 
with 3 j pc 1. As in Fig. 2 all quantities presented 

in this figure correspond to average values over 1 
stellar period (4 th rotation). The Poynting flux 
behaves accordingly, relaxing to L — G.55^ 2 £2 4 /c 3 
for a = 0° (rather than zero), L ~ L0 frft 4 /c 3 for 
a 45° and L — 1 Ap~fl* ft* for a 90° as a — * Q. 
These values are 55%, 65% and 70% of the cor- 
responding FFE values, respectively. For high a 
values (a oc) the corresponding solutions have 


a Poyritiiig flux close to that of the FFE solutions. 
The dissipation power in this case is given by 


E 0 - J J • E dV 

ri<r<r$ 

J E 0 dF 

ri<r<rj 


(n) 


and, as expected, is non-zero even for a - 0, For 
a — Q we get the highest value for Ed t either 
in absolute sense or as a fraction of the Poynt- 
iug flux. These reach 45%, 25% and 12% for 
a — 0 a T a = 45°, a = 90 Q respectively. These are 
higher than the highest values we got using pre- 
scription (B). Their increased magnitude is mostly 
due to the increase in the Ed of t he SFE solutions 
on the current sheet outside the LC. We observe 
also that the dissipation remains at high levels 
even for very high values of the conductivity and 
only the a — 0° case (which is the most dissipative 
one) si lows a clear tendency of decreasing Ed with 
<r . For & > 1 the dissipation energy inside the LC 
decreases significantly while it remains almost con- 
stant outside the LC. We have tried higher resolu- 
tion simulations that allow higher values of & with- 
out any apparent decrease of the dissipation on the 
current sheet This appears to confirm Gruzinov’s 
claim that SFE allows for non-zero dissipation on 
the current sheet even for a -> oo. 

Figure 7 depicts the same quantities as Figs, 3 , 

4 but for prescription (C) Le, SFE and for a c- 8fL 
We observe that the magnetospheric structures in 
this case are very close to those of FFE (especially 
for a > 0°), The magnetic field lines open not 
much farther than the LC and the usual equate- 
rial current sheet develops beyond that, hi the 
third column, which shows the parallel electric 
field, we see that In the effectively time- like re- 
gions near the star the color is “noisy 1 as a result 
of the constant change of its sign during the evo- 
lution, In the next section we are going to dis- 
cuss possible implications of these types of behav- 
ior in these regions. The regions with the "calm" 
(non noisy) electric fields correspond to space-like 
currents. We discern these electric fields mostly 
along the separatrix in the aligned rotator and 
long the undulating equatorial current sheet in 
the non-aligned rotator. We note that for a = 90° 
the entire polar cap region consists of such com- 
ponents. In Fig. S we present the results for the 


same prescription, i.e. (C) - SFE, but for a = 0. 
As we have mentioned in §2, these cases are sig- 
nificantly different from the vacuum solution. The 
magnetic field lines close well outside the LC while 
we are able to discern traces of t he equatorial cur- 
rent sheet, YVe observe again the “noisy 1 ” regions 
near the star (this is the way that SFE deals with 
times like currents, namely as fluctuating space-like 
ones). In the axisymmetric case (a - 0°) we see 
that the value of E\\ as well as the area over which 
it occurs increase for the smaller values of For 
u > the value that reaches near the polar 
cap is not significantly higher that that of Fig. 7, 
However, the new interesting feature in the oblique 
rotators is the development of high E\ \ in specific 
regions of the dosed field lines areas. These con- 
figurations do not imply high dissipation since the 
corresponding currents are very weak there. 

In order to reveal the on-average behavior of the 
effectively time-like regions we considered the av- 
erage values of the fields within one stellar period. 
In Fig. 9 we plot for the SFE (a = 0) the cor- 
responding (average) parallel electric components 
E \\ together with the magnetic field lines (similar 
to the third columns in Figs. 7 and 8). We see that 
the “noisy* behavior has disappeared. As we dis- 
cussed previously one can imagine the existence of 
hybrid solutions (e,g. non lt noisy fl behavior in the 
closed field line regions like that shown in Fig. 9) 
and a “noisy” one in the open field line regions like 
that shown in Fig. 7. 

Finally* we note that for all the solutions pre- 
sented in this section* the regions corresponding 
to high values of J ■ E are traced mostly by the 
high values of J (e.g. near the current layers). 
However, as we observe in Figs. 3-4, 7-8 there are 
regions corresponding to high values of £ !r that do 
not lie near the current sheet. This implies that 
high energy particles supporting the pulsar emis- 
sion may be generated not near the regions of high 
Ed but near those with high values of E \ . 

4, Solutions with J - pc 

As argued above, the time evolution of the 
MHD equations requires a prescription that re- 
lates the current density J and the fields, given 
above by the different forms of Ohm's law. How- 
ever, given that any value of Ej. appropriate for 
pulsars would quickly accelerate any charges to 
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velocities close to c, it would be reasonable to as- 
sume that the relation J ^ pc will be attained 
between the current and charge densities* Moti- 
vated by these considerations, rather than deter- 
mine the current J from its relation to E|j, we 
determine it from its relation to the local charge 
density p. Specifically, we focus our attention to 
solutions with J/pc = L This specific ratio is 
special because it discriminates between magneto- 
spheres consistent with only a single sign charge 
carrier ( J/pc & 1) and those that require carri- 
ers of both signs of charges ( J/pc > 1 ) to ensure 
the counter streaming required in space-like cur- 
rent Rows, Since each how line with J/pc > 1 
requires at some point a source of charges of the 
opposite sign, It is considered that such flows by 
necessity involve the production of pairs and their 
associated cascades. 

One well-studied location for charge produc- 
tion in the pulsar magnetosphere is the polar cap, 
where particle acceleration and pair cascades are 
thought to take place* If charges of either sign 
are freely supplied by the neutron star surface, 
as is thought to be the case in ail but magnetar- 
strength fields (Medm k Lai 2010 ), then a max- 
imum current of J/pc ~ 1 can flow outward 
above the polar cap in Space-Charge Limited Flow 
(SCLF). In this case, an accelerating electric field 
E|l is produced by a charge deficit Ap = (p — pcj) 
that develops because the GJ density is differ- 
ent from the actual charge How along the field 
lines {Arons k Scharlemnmi 1979), This electric 
field then supports a current of primary electrons 
that radiates 7 -rays, initiating pair cascades* The 
pairs efficiently screen the E)] above a pair for- 
mation front (PFF). by downward acceleration of 
positrons* This charge density excess is small rel- 
ative to the total charge* i.e. A p/p = £ <£ 1 
and tiie current that flows into the magnetosphere 
above the PFF has J/pc zz 1 (Harding k Mus- 
limov 2001). Therefore the steady-state SCLF 
acceleration models arc only compatible with a 
small range { J/pc ^ 1 ) of the current to cl large 
density ratio* Beloborodov (2008) demonstrated 
that no discharge occurs if J/pc < I and that the 
E and the discharge generated when J/pc > 1 
or J/pc < 0 (implying that the current direction 
is opposite to the one that corresponds to out- 
ward movement of the Goldreieh- Julian charges) 
are strongly time dependent. This lias been con- 


firmed by models of time-dependent SCLF pair 
cascades (Tiinokhin k Arons 2011 )* 

With this emphasis on the ratio J/pc i, 
we show in Figure 10 the values of J/pc in color 
scale together with the magnetic field lines on the 
poloidal plane for the parameters of the simula- 
tions of prescription (B) given earlier in Figs. 3- 
4 and prescription (C) given in Figs* 7-8. These 
plots can be directly compared to the FFE mag- 
netosphere of Fig. 1 . 

The first two columns of Fig. 10 indicate that 
as a in prescription (B) ranges from large values 
(tj ss 24ft; the FFE has strictly speaking o — oc) 
to 0 (vacuum) the ratio J/pc decreases tending to 
0 for a — ^ 0*° Tins can be seen by inspection of 
Eq, (5), the first term of which corresponds always 
to a time-like current, with a space-like current 
arising only through the effects of the second term. 
Hence there must be some spatial distribution of 
a values below which no space-like region exists 
because the gE\\ term is small enough* 

The next two columns of Fig. 10, showing the 
J/pc ratio of Prescription (C) - SFE, are interest- 
ing because this prescription leads to J/pc 1 = i 
for a — ► 0, while tending near the FFE values 
for o — > oc, but only in the space-like current re- 
gions of the FFE solutions. The time-like current 
regions of the FFE solutions are now only “effec- 
tively time-like”, as discussed by Gruzinov. result- 
ing in the “fluctuating” or “noisy 11 behavior that 
Is apparent there. 

Suspecting that the “fluctuating” (or “noisy 11 ) 
behavior of the effectively time-like current regions 
of SFE is due to the imposed relation between J 
and pc, we modified prescription (B) to search for 
local values of a that would be consistent with 
J/pc= 1 . More specifically, we consider a a vari- 
able in time and space so that the total J provided 
by Eq. (5) is everywhere and at every time step, 
equal to pc. However, we found that there arc re- 
gions where £|j goes to zero for a current J less 
than pc\ In these regions the second term in the 
expression for ./. namely <riE?j|, tends to the second 
term of the FFE current expression (Eq. 3) be- 
ing unable to comply with the J pc constraint. 
These regions demand extremely high values of 
a in order to achieve the J/pc 1 relation due 


f ' The oblique rotators have J/^c ratios near 0 for much lower 
values of & than those of the first column of Fig. 10* 
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to the vanishing values. We have dealt with 
this numerically by setting a maximum value of a 
(-- 50Q) that provided a realistic suppression of 
the E\\ in these regions. 

Figure 11 shows the results of such an approach. 
The first three columns provide respectively the 
values of the poloidal current J p , charge density p 
and parallel electric component E [j and are simi- 
lar to those presented in Figs. 3, 4> 7* 8, while the 
fourth column shows the values of the Jj pc ratio 
in color scale together with the poloidal magnetic 
field lines. The corresponding L and E& values 
as well as their ratios Eo/L have been denoted 
in Fig. 2 by the horizontal dashed lines. Not sur- 
prisingly, these solutions are similar to those of 
prescription (C) - SFE in the a — > 0 regime (see 
Fig. 8). The essential difference between these 
two approaches is that, in distinction with SFE, 
this approach can handle time-like currents with- 
out the need - like SFE - to do so only on average 
(see blue and purple color regions in the fourth 
column). 

To allow our solutions even larger flexibility in 
attaining J = pc we have looked for prescriptions 
that divorce J from Ey, replacing the second term 
of Eq. (5). with a term of the form 

(12) 

where / is a scalar quantity having a local value 
that ensures the local value of the ratio J/pc to 
be equal to 1. The difference from the previous 
approach and that of prescription (B) is that this 
term is not related anymore to the value of E\\. 
Other than that it is directed along (or against for 
/ < 0) the direction of the magnetic field. The 
sign of / is determined so that the direction of J ( [ 
coincides with that of E||. The results of these sim- 
ulations are shown in Fig. 12 the panels of which 
are similar to those of Fig. 11. The global struc- 
ture in this case is very similar to that of the simu- 
lation presented in Fig. 11 except for the time-like 
regions (especially those near the star) that now 
exhibit the "noisy" behavior of the SFE shnula- 
t i o ns . Howe ver . the ratio J/ pc is \ 10 w every w b ere 
equal to 1. Our solution has therefore the required 
property (J/pc = 1) at the expense of producing 
a “noisy" behavior indicative of time variations at 
these limited spatial scales. 

This oscillatory behavior (both here and in the 
SFE cases) is a rather generic feature that occurs 


in the time-like current regions for values of J/pc 
larger than those needed to “kill” the local Ey im- 
plied by the field evolution equations; both within 
the prescription of Eq. (12) and within SFE, at 
these higher values, ./ and Eu end up pointing in 
different directions and an oscillatory behavior en- 
sues. Smaller values of J/pc , cannot drive Ej| to 
zero and allow the magnetosphere to attain steady 
state, leaving a uon-oscillatory component of Eg 
along the direction of J|p Viewed this way, the 
idea) condition E - B 0 is one that requires a lo- 
cal delicate balance between its parameters (e.g. 
J,p). If, for whatever reason or necessity, this bal- 
ance breaks, the magnetosphere develops locally 
an oscillatory behavior. The robustness and the 
true evolution of this behavior will be determined 
by the microphysical plasma properties as they 
mesh within the global framework of the magne- 
tosphere electrodynamics. 

The conclusion of this exercise is that solutions 
with a specific value of J f for example J — pc ev- 
erywhere, can be implemented with a sufficiently 
general expression for J y , provided that one is will- 
ing to accept locally time varying solutions. 

Note that in this prescription, as well as in 
SFE, the current is forced to follow the relation 
J = pc or J > pc (for a ^ 0 in SFE) even in the 
closed field line regions producing even there this 
“noisy” (fluctuating) behavior. Although apply- 
ing this J — p prescription there is not forbidden a 
priori, it would be reasonable to consider that this 
condition is actually not valid in the closed field 
lines and that a different condition should be used 
there. However, the implementation of a different 
prescription for these regions remains cumbersome 
because it demands that we determine whether the 
magnetic field lines at each grid point are actually 
open or closed and then apply for them this dif- 
ferent condition. Nevertheless, even if we were to 
enforce a prescription that removes the fluctuat- 
ing i noisy) behavior in the closed field line regions, 
this should not affect the similar behavior on the 
open field line regions. 

5, Discussion and Conclusions 

Up to now, the pulsar radiation problem 
has been studied considering either the vacuum 
(Deutsch 1955) or the FFE solutions (Contopou- 
los et ah 1999: Spitkovsky 2006; KaJapotharakos 
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& Coutopoulos 2009) for the structure of the un- 
derlying magnetosphere. The vacuum solutions 
have analytical expressions while the FFE solu- 
tions have only been studied numerically in the 
last decade. However, neither solution is com- 
patible with the observed radiation; the former 
because it is devoid of particles and the latter be- 
cause it precludes the presence of any Ej* The 
physically acceptable solutions have to lie some- 
where in the middle. As of now, there are no 
known solutions that incorporate self-consistently 
the global magnetospheric structure along with 
the microphysics of particle acceleration and ra- 
diation emission. The goal of the present work 
has been to explore the properties of pulsar mag- 
netospheres under a variety of prescriptions for 
the macroscopic properties of the underlying mi- 
crophysics in anticipation of a more detailed fu- 
ture treatment of these processes. It should be 
noted that these processes lie outside the purview 
of the equations used to evolve the electric and 
magnetic fields. Nonetheless, for most models 
examined the dissipation power does not exceed 
10-20% of the spin-down one; this value is con- 
sistent with the observed radiative efficiencies of 
the millisecond pulsars f Abdo et ah 2009a) t that 
are the ones reaching the highest efficiencies of 
all pulsars, suggesting that energetically, realistic 
magnetospheres may be well approximated by our 
models. 

In order to explore the non-ideal E B ^ 0 
regime we considered several prescriptions dis- 
cussed in §3, 4 and analyzed their properties in 
detail (Figures 2-9)- Despite their intrinsic in- 
terest, the solutions discussed above still need to 
be examined more closely for consistency with the 
microphysics responsible for the set up and clo- 
sure of the corresponding magnetospheric circuits, 
which arc characterized by a considerable range in 
the distribution of the ratio Jj pc. On the other 
hand, there axe studies arguing that it is essential 
that the ratio J / pc be near 1 in order for steady- 
state pair cascades to be generated near the po- 
lar cap. Besides that, solutions corresponding to 
values J/pc < I can be considered as the range of 
con figurations without pair production, since such 
solutions can be supported by charge separated 
flows. 

With the value J/pc 1 given this particular 
significance, we have pursued a search for solu- 


tions that are driven by adherence to this condi- 
tion in as broad a spatial range as possible. As 
noted in the previous section we did so by search- 
ing for solutions of variable conductivity u , or by 
employing expressions for the parallel current un- 
constrained by the values of the parallel electric 
field. We found that arbitrary values of <r, while 
allowing J/pc = 1 in most areas, in certain others 
can drive E|j to zero sufficiently fast to restrict J 
to values J < pc. As rioted in the previous sec- 
tion, a more general expression for J not related 
to E|| can produce J/pc = l everywhere but only 
at the expense of fluctuations in regions where the 
value of J/pc is larger than that necessary to null 
the local E;| implied by the field evolution equa- 
tions, These results indicate the importance of 
the value of the ratio J/pc vis-a-vis the presence 
of pair cascades and steady-state emission in pul- 
sar magnetospheres, 

The actual behavior of the magnetosphere is 
likely to be more complicated than that described 
by the current density prescriptions provided in 
this work (Eqs, 5, 6 and 12). Physics beyond those 
of the magnetic and electric field evolution equa- 
tions will introduce novel time scales that will de- 
termine the response of the current density to the 
parallel electric field component and vice versa. 
The implications of such an oscillatory behavior 
which involves an electric field that alternates di- 
rection t may in fact provide an account of the pul- 
sar coherent radio emission, as such a behavior is 
the essential ingredient of cyclotron maser emis- 
sion (Melrose 1978; Levinson et ah 2005; Melrose 
2006; Luo Melrose 2008), To the best of our 
knowledge, the above argument Is the first ever to 
suggest a relation between the global properties 
of pulsar magnetospheres and their coherent radio 
emission. 

These more general models of non- ideal magne- 
tospheres discussed above are possibly related to 
the behavior and properties of intermittent pul- 
sars, which may, if more of them are discovered, 
provide meaningful constraints on their radiation 
mechanism and magnetospheric dynamics. Inter- 
mit tent pulsars are pulsars found in "on" and 
"ofT radio states (Lyne 2009). The interesting 
aspect of these objects is that they exhibit dif- 
ferent spin-down rates in their "on" and "off* 
states; even more interesting is the fact that they 
apparently "remember' the spin-down rate they 
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had in the previous cycle, as they return to the 
same spin-down rate when switching from “on" 
to “oflF” and vice versa. The ratios of the u on n 
to “ofT spin-down rates are 1.5 -1,7 for pul- 
sars PSR B 1931 +24 (Kramer et a 1. 2006) and 
PSR J 1832+0029 (Lorimer et aL 2006: Lyne 2009), 
while it is 2.5 for the more recently observed 
PSR J 1841-0500 (Camilo et &L 2011). Interest- 
ingly these values are well below the minimum ra- 
tio 3 between the FFE and the vacuum spin-down 
rates. Based on tills observation, Li et ah (2011) 
argued that the “ofP state can be the vacuum 
solution while the “on” state simply a resistive so- 
lution corresponding to a certain value of cr. 

An alternative view is to Identify the “on” state 
with magnetosperes with J/pc > 1 and the “off” 
one with those of J/pc < 1, given the necessity 
of the presence of pairs (accompanied by the ra- 
dio emission) in the former and their likely ab- 
sence (and also of radiation) in the latter. In fact, 
the ratio of spin-down rates between the FFE and 
the J j pc — 1“ magnetospheres is within 1.4- 1.8 
(see Figs, 2 and 6), with the higher values corre- 
sponding to the lower values of a. This would be 
consistent with the observations of the first two in- 
termittent pulsars. However, this is not consistent 
with the spin-down ratio of the third intermittent 
pulsar whose “on” and “off” states may be magne- 
tospheres with different values of o\ The discovery 
of a larger sample of intermittent pulsars (espe- 
cially if one could also put some limits on their 
inclination angles) would be extremely useful in 
constraining magnet os pheric models, For exam- 
ple if the ratio of the a on/ofT spin-down rate was 
found not to exceed the value 3, one would have to 
exclude the vacuum solution being the “off' state, 
since for sufficiently small a this ratio becomes ar- 
bitrarily large, 

The solutions we presented in this work, besides 
their intrinsic interest, can serve to produce model 
pulsar "—ray light curves in a way similar to that 
of Contopoulos & Kalapotharakos (2010): Bai & 
Spitkovsky (2010a); Venter et u3. (2009): Harding 
et ah (2011). Such models will provide a means of 
connecting the structure and physics of magneto- 
spheres to observations, since they constrain the 
pulsar magnet os pheric geometry as well as the re- 
gions of anticipated particle acceleration and even 
the values of the accelerating electric fields. These 
will therefore provide independent tests of these 


models for real pulsars. 

The next level of this study should include mi- 
crophysics at a level sufficient to allow some feed- 
back on the global solutions i,e. the production 
of pairs whose effects are included in the com- 
putation of the local E\\ self-consistently. It will 
be of interest to examine whether such solutions 
are consistent with the charge separated flows or 
whether they require an inherently time dependent 
magnetosphere, a feature that could be tested ob- 
servation ally. The time dependence should lead to 
oscillations of the pulsar emission with frequency 
of order v Q/2n for emission near the LC and 
v :§> 0/2tt for emission near the polar cap. We 
have searched and found not these oscillations in 
the X-ray domain of several pulsar light curves. 
Either they were not present at the objects we 
searched or they lose their coherence at these much 
lower energies. Nonetheless, we intend to study 
the possibility of the presence of such oscillation 
in a future work. 

We would like to thank the anonymous referee 
for constructive comments that helped us to im- 
prove this paper. 
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Fig, 1. — The FFE solutions for the inclination angles a indicated in the figure. Each panel shows the 
poloidal magnetic field lines (white colored lines) arid the ratio Jj pc in the color scale shown in the figure. 
We note that the color scale saturates below the value 0.5 and above 1.5 so that details around the value 
I become evident. The dashed vertical lines denote the Light Cylinder (LC). The length unit is equal to 
Rlc* We observe that the magnetic field lines open beyond the LC, The fraction of the polar cap region 
with space-like currents (J/pc > 1) increases with a so that, while in the aligned rotator almost the entire 
polar cap is filled with time- like currents ( J/pc < 1), in the perpendicular rotator the entire polar cap is 
filled wit! i space-like currents (J/pc > 1). 
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Fig, 2* (a) The Poynting fluxes L (measured on 

the surface of the star) as a function of a in log- 
linear scale for prescription (B), Red, green and 
blue colors correspond to a - 0° (aligned), a = 45° 
and a — 90° perpendicular rotator, respectively. 
The horizontal solid line elements denote the L 
values corresponding to the vacuum (lower value) 
and the FFE (higher value) solutions* Note that 
the L value for a 0 s is similar to the vacuum one 
for rr < 0.3ft (ft is the angular frequency of the 
star) and it reaches that of the FFE solution only 
for much higher a values: (b) The dissipation en- 
ergy rate Eq integrated over the volume bounded 
by radii r t = 0.3 Rlc and r '2 — 2.51?lc as 
a function of <r: (c) The fraction Epf L never ex- 
ceeds the values 10%-209F while for a — > ► 0 and 
<t — y x; it goes towards 0* The dashed horizontal 
lines in ail three panels denote the values corre- 
sponding to simulations with J/pc - 1 (see §4, 
5}. 
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Fig. 3. Tlie non- Ideal magnetospheric solutions in the poloidal plane (ft. SI) for prescription (B) and tor 
a high a value (a = 2417). Each row shows the solutions corresponding to the indicated inclination angles 
a. The left-hand column shows the poloidal current modulus (in color scale) together with the streamlines 
of the poloidal current. The middle column shows the charge density (in color scale) together with the 
field lines of the poloidal magnetic field. The color ranges purple-green and green-red indicate negative and 
positive charge density, respectively. The right-hand column shows the parallel electric field component E 
(in color scale) together with the lines of the poloidal magnetic field. The color ranges purple-green and 
green-red indicate antiparallel and parallel directions of E \ (relative to the magnetic field) respectively. Note 
that the color representation for E saturates beyond the absolute value 0,2 j Blc/c where Bu: is the value of 
the stellar magnetic dipole field at the distance /?lc- The structure shown in the first two columns is quite 
similar to that of the FFE solutions ( the main difference from the FFE solution can be seen in the third 
column). 
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Fig. 4. — Similar to Fig, 3 but for a much lower a value a = L5f2), The main magnetospheri c features such 
as separatrices and current sheets can still be observed, even though they now appear much weaker. The 
magnetic field lines close well beyond the LC and the parallel electric field components reach higher values 
than those of Fig. 3 
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Fig* 5. - Similar to the third columns of Figs. 3* 
4 but for the prescription (A) and for the h — 0.75 
value. For this value of h the corresponding L 
values are the same to those of t he simulations 
presented in Fig. 4. The structure of these solu- 
tions seems very similar qualitatively to those of 
Fig. 4. 
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Fig. 6* — Similar to Fig* 2 but for prescription 
(C) - SFE* Tills prescription does not cover the 
entire spectrum of solutions between the vacuum 
and FFE ones since for rr — > 0 we still get a dis- 
sipative configuration rather than the vacuum so- 
lution* However* as a x* the corresponding 
Poyuting flux values tend to those of the FFE so- 
lutions. The dissipative energy loss rate Eq ex- 
hibits a maximum value for <7 — 0 and remains 
non-zero even for the high a values due to the per- 
sistent dissipation on the current layer outside the 
LC. This result seems to confirm Gr uzin ov who 
argued that the SFE solutions for or ■=* x present 
finite dissipation on the current layer* 
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Fig* 7. — - Similar to Figs* 2, 3 but for prescription (C) - SFE and for a high u value (a = 8Q)* The Poynting 
[lux of these solutions tend to those of the FFE ones. Nevertheless, we still see magnetic field lines closing 
well outside the LC (especially for a 0°). The electric current in the SFE prescription is everywhere space* 
like i JJ pc >1). However, there are regions that are effectively time- like (as defined in the text These are 
the noisy regions in the polar cap vicinity (for a / 90 : ; seen in the third column* Their noisy character is 
due to the continuously changing direction of the electric field* 
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a= 90° 


Fig, 8,— Similar to Fig, 7 but for a = 0 value. The modulus of the parallel electric field components 
increases along with the area of their influence. The closed field line regions of the oblique rotators show 
significant parallel electric field components: however, the dissipation there is insignificant because the value 
of the product J ■ E is small. 
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Fig. 9, Tlie third column of Fig. S considering 
the average values of the fields within a stellar pe- 
riod. The " noisy" behavior shown in Fig. 8 is now 
disappeared. These regions reveal their effectively 
time-like behavior* 
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Fig. 10. — The ratio J/pc (in color scale) together with the poloidal magnetic field lines for the solutions 
presented in Figs. 3, 4, 7, 8 as indicated in the Figure. The ratio J/pc in prescription (B) decreases in 
general as a decreases and goes towards 0 for a -*> 0. In prescription (C) - SFE the ratio J/pc goes to 1 
(null current) as <r — * 0 while it is higher than 1 for t? > 0. However, in both these cases there are regions 
where the current has on average time-like (J/pc < 1) behavior (effectively time-tike: see Fig, 9). 
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Fig, 11,— The solutions employing prescription (B) with the <r — values adapted locally to achieve J/pc as 
close to 1 as possible. The first three columns are similar to those of Figs. 3, 4, 7, 8 while the fourth column 
shows the corresponding ratio J/pc in color scale. The main features (e.g. current sheets, separatrices) of 
the FFE solutions are still discernible. We observe parallel electric field components across the polar caps 
and along the separatrices. For the oblique rotators parallel electric field components exist also in parts 
of the closed field line regions. However, the latter regions show low dissipation energy rate Eq since the 
corresponding product J E is small* In the fourth column we observe regions with J/pc < L In these 
regions a current density J smaller than the value pc is enough so that the corresponding parallel electric 
field component vanishes. This kind of solutions are the upper limit of the solutions that can be considered 
that are supported by charge separated flows. 
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= 90° 


Fig, 12.” Similar to Fig. II but with the parallel current density independent of the modulus of E (see 
Eq, 12). As shown in the fourth column, tins prescription achieves Jf pc = 1 in the entire magnetosphere. 
This leads to an oscillating behavior for E in the time-like current regions of the simulation presented in 
Fig. II* similar to that of SFE {Fig. 8). We note that this effect is more general and appears wherever 
the ratio Jf pc is enforced to have values above those that are sufficient to cancel out the parallel electric 
component. This behavior is the result of the mutual response between the electric held and the current 
density and may be crucial for the coherent radio emission. 
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